Astronomy & Astrophysics manuscript no. ms 


©ESO 2013 


January 20, 2013 





Secular spin-down of the AMP XTE J1 751-305 

A. Riggio 1 - 2 , L. Burderi 2 , T. Di Salvo 3 , A. Papitto 1 - 2 , A. D'Ai 3 , R. Iaria 3 and M. T. Menna 4 

1 INAF/Osservatorio Astronomico di Cagliari, localita Poggio dei Pini, strada 54, 09012 Capoterra, Italy; e-mail: 
ariggio@oa-cagliari . inaf . it 

2 Universita di Cagliari, Dipartimento di Fisica, SP Monserrato-Sestu km 0,7, 09042 Monserrato (CA), Italy 

3 Dipartimento di Scienze Fisiche e Astronomiche, Universita di Palermo, Via Archirafi 36, Palermo, 90123 Italy 

4 Osservatorio Astronomico di Roma, Sede di Monteporzio Catone, Via Frascati 33, Roma, 00040 Italy 



o 



x 

6 

(N 
> 

oo 



in 
o 



X 



Abstract 

Context. Of the 13 known accreting millisecond pulsars, only a few showed more than one outburst during the RXTE era. XTE J175 1- 
305 showed, after the main outburst in 2002, other three dim outbursts. We report on the timing analysis of the latest one, occurred on 
October 8, 2009 and serendipitously observed from its very beginning by RXTE. 

Aims. The detection of the pulsation during more than one outburst permits to obtain a better constraint of the orbital parameters and 
their evolution as well as to track the secular spin frequency evolution of the source. 

Methods. Using the RXTE data of the last outburst of the AMP XTE J175 1-305, we performed a timing analysis to improve the 
orbital parameters. Because of the low statistics, we used an epoch folding search technique on the whole data set to improve the local 
estimate of the time of ascending node passage. 

Results. Using this new orbital solution we epoch folded data obtaining three pulse phase delays on a time span of 1.2 days, that we 
fitted using a constant spin frequency model. Comparing this barycentric spin frequency with that of the 2002 outburst, we obtained a 
secular spin frequency derivative of -0.55(12) x 10~ 14 Hz s -1 . In the hypothesis that the secular spin-down is due to a rotating magneto 
dipole emission, consistently with what is assumed for radio pulsars, we estimate the pulsar's magnetic dipole value. We derive an 
estimate of the magnetic field strength at the polar cap of B P c = 4.0(4) x 10 8 Gauss, for a neutron star mass of 1.4M , assuming the 
Friedman Pandharipande Skyrme equation of state. 

Key words, stars: neutron - stars: magnetic fields - pulsars: general - pulsars: individual: XTE J175 1-305 - X-ray: binaries. 



1. Introduction 

XTE J175 1-305 is one of the accretion powered millisec- 
ond X-ray pulsars (AMPs) that showed more than one out- 
burst in the RXTE era. Re current outbursts were also ob- 
served in SAX J1808 4-36 58 (Idi Salvo et al.l2008l:lBurderi et al 
2009t lHartman et al.ll2009h IGR J00291+5934 dGallowav et al 



2005 



120081: lPatrunoll2010t lHartman et ail 1201 U IPapitto et al 



20101) and recently in N GC 6440 X-2 (lAltamirano et al.ll2oTol) 



and Swift J1756.9-2508 dPatruno et al.ll2010h . XTE J175 1-305 
was detected for the first t ime by RXTE on April 3, 2002 
(iMarkwardt et al.ll2002llM02l henceforth). This outburst was the 
brightest and longe st of the four showe d by XTE J1751-305, 
permitting !m02| and lPapitto et al.1 d2008l) to obtain a full orbital 
solution. 

The second o utburst, was spotted by INTEGRAL 
dGrebenev et al.1 120051) on March 28, 2005 and lasted ~2 
days with a peak flux which was ~ 14% of the one reached 
during the first outburst. Unfortunately the RXTE PCA in- 
strument was not in the p roper mode to d etect pulsations 
during the first observation (ISwank et al.ll2005l) . The follow-up 
observations was done in event and single bit modes, but the 
source already fainted below detection. The attribution of this 
outburst to XTE J1751-305 is not certain, since pulsations were 
not detected and the INTEGRAL IBIS/ISGRI source position 
is compatible with the position of at least one other s ource, 
IGR J1751 1-3057 (IPapitto et alJ20!otlRlggio et alJl20Tlh . 

The third outburst was detected on April 5, 2007 by 
RXTE with a peak flux of 18% of the first outburst 
(IMarkwardt & Swankl2607tlFalanga et al.l2007l) . very similar to 



the second outburst. In the subsequent pointed observation by 
RXTE the source became too dim to detect pulsations. In this 
case the source id entification is certain, th anks to a simultaneous 
Swift observation (Markwardt et al. 20Q7|). 

The latest outbu rst was first spotted with INTEGRAL 
dChenevez et al. 2009). Fortunately, it occurred while RXTE was 
observing the last di scovered AMP IGR J 17 5 1 1-3057, very near 
to XTE J1751-305 (IMarkwardt et al.l2009l) . The detection of the 
435 Hz X-ray pulsation and a following Swift observation con- 
firmed that the source in outburst was XTE J1751-305 and not a 
re-brightening of IGR J1751 1-3057. In this paper we report on 
the timing analysis of this outburst. 



2. Observation and Data Analysis 

In this work we analyse RXTE PCA observations of XTE J175 1- 
305. We use data f rom the PCA (proportional counter array, see 
iJahoda et al J 12006) instrument on board of the RXTE satellite 
(Obsld 94041 and 94042). We used data collected in event pack- 
ing mode, with time and energy resolution of 122 fj.s and 64 en- 
ergy channels respectively. Although the XTE J1751-305 pulsa- 
tion was detecte d only during observati ons performed on 8 and 
9 October 2009 (IMarkwardt et al.ll2009l) . the data analysed here 
cover the time span from 6 October 2009 to 22 October 2009, 
as shown in Figure Q] This allows to precisely determine any 
contribution to the observed emission attributable to the AMP 
IGR J1751 1-3057, that was still in outburst during the RXTE 
observation, and to the galactic ridge. We obtained the energy 
band which maximises the signal to noise ratio comparing the 



1 



A. Riggio et al.: Secular spin-down of the AMP XTE J1751-305 



45 



40 \ . 




1 1 1 1 1 ^ 

-5 5 10 15 

Time (MJD - 55112.0) 

Figure 1. PCU 2 count rate (2-18 keV), subtracted of its back- 
ground, is reported as a function of time in the period from 
2 October 2009 and 22 October 2009. During the first days, 
the last phase of the IGR J1751 1-3057 outburst is visible. The 
flux re-brightening is caused by the onset of the XTE J 1751- 
305 outburst, which lasted less than two days. The remaining 
days show the constant flux due to the galactic ridge. The su- 
perimposed model is the best-fit using a piecewise linear func- 
tion. Since we are interested in determine the background due to 
IGR J1751 1-3057 and the galactic ridge, we excluded from the 
fit the XTE J1751-305 outburst. See the text for more details. 



source X-ray spectrum and the background. We chose the 2-18 
keV energy band. We corrected the photon arrival times for the 
motion of the earth-spacecraft system with respect to the solar 
system barycentre and reported them to barycentric dynamical 
times at the solar system barycentre using the faxbary tool (DE- 
405 solar system ephemeris), adopting the Chandra source po- 
sition rep orted by lM02i The uncertainty on the source position 
quoted by M02j is 0.6", 90% confidence level, as shown in Table 

m 

We obtain a first estimate of the mean spin frequency con- 
structing a Fourier power density spectrum of the 3.2 ksec expo- 
sure ObsID 94041-01-04-08 data calculating 53 power spectra 
from 64-s long data segments (2~ n bin size) which were aver- 
aged i nto one power spectrum. As reported in iMarkwardt et al.l 
(2009), we found a signal at ~ 435.32 Hz. No conclusive evi- 
dence of pulsations in the following observations was found in 
this preliminary step. 

2.1. Determination of the local T* 

A timing analysis to obtain an orbital solution at the time of the 
latest outburst from XTE J175 1-305 was not possible because 
of the weakness of the pulsation. However, it is still possible to 
correct the time series for delays induced by the orbital motion 
adopting t he orbital paramete rs of th e April 2002 outburst esti- 
mated by iPapitto et al.l (120081) (|P08l hereafter). Propagating the 
uncertainty in the orbital period P 01 -b given by |P08l along the ~ 7 
years separating the 2009 outburst from the 2002 one, resulted 
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Figure 2. y 2 maxima obtained from an epoch folding search on 
the 2009 data corrected for the orbital modulation, varying each 
time the epoch of passage through the ascending node with 1 s 
step. A total of 2546 orbital solutions were tried, exploring all 
the possible values for T*. 



in an uncertainty on the time of passage through the ascending 
node, T*, in the 2009 outburst of ~ 186 s at lcr confidence level, 
about ~ 7% of the orbital period. Moreover the presence of an 
orbital period derivative P 01 -b might introduce a further shift on 
T*. The true local value of T* can thus be significantly differ- 
ent from the nominal value obtained by propagating the orbital 
solution provided by P08. 

To determine the best local orbital solution, we make the rea- 
sonable hypothesis that the best local set of orbital parameters is 
the one which gives the best signal to noise ratio, that is, in our 
case, the highest y 2 value in an epoch folding search (see e.g. 
iKirsch et alJ l2004). We restrict our search to just one orbital pa- 
rameter, T*, which is the parameter with the largest uncertainty 
during the 2009 outburst. To explore all the possible values for 
T*, the orbital period being 2546 s, we produced 2546 differ- 
ent time series from the data of the 2009 outburst, which were 
corrected for the orbital modulation. For each of these time se- 
ries the adopted orbital parameters were the same, except for 
T*, which is varied in steps of 1 s. We then performed an epoch 
folding search for the spin period on each of the 2546 time series 
using 32 phase bins to sample the signal. In Fig. [2] we show the 
maximum of the x 1 obtained from the epoch folding search of 
each time series as a function of the T* adopted to produce the 
time series on which the folding search has been performed. We 
fitted the^ 2 maxima curve with a model consisting of a constant 
plus a Gaussian. A clear peak at AT* -1 10 s is evident, well 
with in the lcr confidence level of the orbital solution given by 
P08 propagated to the 2009 outburst. Thus we derived the time 
of passage from the ascending node during the 2009 outburst as 
T* = T* + AT* 

Adopting the new value of T*, we barycentered all the data 
covering the 2009 outburst of XTE J1751-305 and performed 
again an epoch folding search of the spin period. The improve- 
ment of the orbital solution allowed us to detect the pulsation in 
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Table 1. Orbital and Spin Parameters for XTE J1751-305. 




Parameter 



Value 



-125 



-115 -110 -105 -100 
T* - T* P08 (seconds) 



Figure 3. Maximum of the x 2 obtained in an epoch folding 
search around the expected spin period as a function of the T* 
adopted to correct the time series for the delays induced by the 
orbital motion. In this figure we show the result of the final 
search performed on all the three observations for which the pul- 
sation was detected after having corrected for the orbital motion 
using the value T* as estimated from Fig. [2] A T* step of 0.15 s 
was adopted and 200 T* were tried. The best-fit model, constant 
+ Gaussian, is also shown (dashed curve). 



two other observations, corresponding to ObsID 94041-01-04- 
04 (MJD 55112.335, with an exposure of 1.3 ksec), and ObsID 
94042-0 1 -02-00 (MJD 55 1 1 3 . 1 65, with an exposure of 3 . 1 ksec, 
see Tab.[2]for details). 

We applied the procedure described above considering all the 
three datasets in order to further improve our new T ( * 9 measure. 
We adopted a finer step in T* of 0. 15 s, covering an interval of 30 
s around the new value of T*. Again, we fitted the^ 2 maxima vs. 
T* curve with a model consisting of a constant plus a Gaussian, 
as it is shown in Fig. [3] In this way we were able to obtain a 
precise measure of T*; the final value for T* is reported in Table 
Q] The evaluation of the uncertainty on T* is discussed in the 
next section. 



2.2. Error estimates on T* using Monte Carlo simulations 

The folding search technique described above and adopted to ob- 
tain our refined measure of T* does not provide a straightforward 
determination of the uncertainty on this parameter. To overcome 
this problem we performed a Monte-Carlo simulation. We gen- 
erated 100 datasets with the same exposure, count rate, pulse 
modulation and orbital modulation observed in the real data. For 
each of these dataset we applied the same procedure as described 
in the previous section to obtain a measure of T* . The confidence 
level at 68% (lcr) for the T* parameter corresponds to 1.05 s, 
which is about 30 times the corresponding confidence interval 
obtained by a fit with a gaussian of the x 2 maxima vs T* curve. 
We therefore adopt 1.05 s as our best estimate of the 1 <x uncer- 
tainty on the T* measure. 



RA (J2000) 
Dec (J2000) 

Projected semi-major axis a x sin i (lt-ms) 
Ascending node passage, T* (MJD) at 

2002 outburst 
2009 outburst 

Orbital period, P mh (s) 
Eccentricity, e 
Reference epoch, To (MJD) 
Mean spin frequency, v (Hz) 
Secular spin frequency derivative, 

v sec (Hz s~') 



17 h 51 m 13?49(5) 
-30°37'23':4(6) 
10.125(5) 

52368.0129023(4) 
55111.000647(12) 
2545.342(2) 
< 1.3 x 10- 3 
55112.0 
435.31799256(23) 

-0.55(12) x 10-' 4 



References. " iMarkwardt et"aT1(l2002h : b iPapitto et aT] d2008h 

Notes. Errors are at lcr confidence level, upper limits are given at 95% 
confidence level. Errors on source position are given at 90% confidence 
level. Times are referred to the barycentre of the Solar System (TDB). 



2.3. Timing analysis 

We barycentered our data of XTE J1751-305 using our refined 
orbital solution, and performed an epoch folding search to de- 
termine a mean spin frequency for the 2009 outburst. From the 
best-fit of the^f 2 curve we obtained a value of 435.31799237(19) 
Hz. We evaluated the frequency uncertainty using the Monte- 
Carlo simulated data described above. We found that the error 
determined in this way was an order of magnitude bigger with 
respect to the error determined fitting with a gaussian the cen- 
troid position of the epoch folding search curve. 

Adopting this mean spin frequency value, we epoch folded 
the three observations during which pulsations could be de- 
tected over 1000 s long intervals, considering 16 bins to sam- 
ple the signal (see Table|2]for details). In this way we obtained 
7 folded pulse profiles. We performed an harmonic decomposi- 
tion of each pulse profile. The fundamental and the first over- 
tone were significantly detected. The fundamental was signifi- 
cantly detected on 6 folded pulse profiles while the first overtone 
only in one profile. We fitted the pulse phase delays with a con- 
stant plus a linear term, representing a constant (mean) spin fre- 
quency model. From the fit we obtained a mean spin frequency 
of 435.31799256(22), with a final x 2 /dof of 1.55(4 dof). This 
value is in perfect agreement with the value obtained with the 
epoch folding search. Note that the uncertainty on the pulse fre- 
quency obtained from the fitting of the pulse phase delays fit is 
nearly equal to the uncertainty on the pulse frequency obtained 
with a folding search estimated with the Monte Carlo simula- 
tions, so that both procedures give consistent results. The pulse 
phase delays and the best fit line are shown in Figure |4] The 
pulse profile obtained by folding all the data is shown in Figure 


To correctly determine the fractional amplitude, an estimate 
of the background and the contribution from sources contam- 
inating the field of view is mandatory. In this case, the major 
contribution to the background is due to the emission of the 
AMP IGR J1751 1-3057 and of the galactic ridge. While the 
roughly constant contribution from the galactic ridge is of ~ 7.5 
cts s _1 in the considered energy band (Pap itto et al.ll2.Ql6T) , the 
IGR J175 1 1-3057 contribution can only be estimated extrapolat- 
ing its flux decay trend just before the XTE J1751-305 outburst 
onset. We fitted the 2-18 keV IGR J1751 1-3057 X-ray light 
curve with a linear model in the time interval from 55108.92 to 
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Table 2. Analysed Obs Id exposures and fractional amplitudes. 



0.4 0.6 0.8 1 
Time (MJD - 55112.0) 



1.4 



Figure 4. Pulse phase delays of the fundamental for the three 
observations in which the pulsation is detected. Each phase point 
is obtained folding on ~ 1000-s long time intervals and using a 
spin frequency v = 435.31799237, which is the value obtained 
with the epoch folding search technique. The best-fit constant 
spin frequency model is also shown (dashed line). 



551 1 1.46 MJD, as it is possible to see in FigureQ] In Table|2]we 
report the obtained fractional amplitudes for each of the three ob- 
servation. These values are corrected with respect to the instru- 
mental background as well as the estimated contribution from 
IGR J1751 1-3057 and the galactic ridge, but are still strongly 
dependent from the model adopted to describe the IGR J 175 1 1— 
3057 X-ray light curve. 

The uncertainty on the position of the source quoted b y M02 
is of 0.6"(90% confidence level, see lAldcroft et al.ll2000l) . while 
0.37" is the confidence interval corresponding to 1 oQ Because 
of this, a systematic uncertainty arises on the spin frequency ob- 
tained by fitting the pulse phase delays. These systematic un- 
certainties are cr v S y S = 4.8 x 10~ 7 v 3 cr" (1 + sin 2 /3) 1/2 Hz 
(Burd eri et al.l l2007). where V3 is the spin frequency in units of 
1000 Hz, cr" is the positional error circle in units of arcsec, and 
f3 refers to the ecliptic latitude of the source (for XTE J 175 1-305 
A = 268.097281° and j3 = -7.198364°). With these values we 



have cr. 



v sys 



: 7.8 x 10~ 8 Hz. 



Combining in quadrature this systematic error with the sta- 
tistical error of 2.2 x 10~ 7 Hz we find a final error on the spin fre- 
quency of o-y = 2.3 x 10 1 Hz. Thus the value of the average spin 
frequency during the 2009 outburst is v 09 = 435.31799256(23) 
Hz. 



3. Discussion 

We obtained a refined orbital solution and a precise estimate of 
the spin frequency of the AMP XTE J1751-305 from a timing 
analysis of the RXTE data during its 2009 outburst. 



Obs ID 


Start time 


Exposure 


fractional amplitude 




(MJD) 


(ks) 


(%) 


94041-01-04-08 


55112.052 


3.2 


7.9+ 1.0 


94041-01-04-04 


55112.335 


1.3 


7.6 ± 2.2 


94042-01-02-00 


55113.165 


3.1 


7.1 ± 2.0 



Notes. For each Obs Id in which pulsations at the frequency of 
XTE J175 1-305 have been detected in this work, the start time, the 
exposure and the fractional amplitude (corrected for the contribution 
of IGR J1751 1-3057, galactic ridge, and instrumental background) are 
reported. 



3.1. Orbital period evolution 

Although a measure of P 01 -b is impossible with only two mea- 
surements of T*, we can derive an upper and l ower limit to the 
P 01 b using the full orbital solution given bv lP08l for the 2002 out- 
burst an d our measure of T* for the 2009 outburst. Using eq. 1 
given in Burd eri et al.l (j2009), we can solve that expression to 
derive P or b 



Pnrh = 



ATP, 



orb 



AT* (AO 
N 



-AP n 



(1) 



where P or b is the orbital period measured by P08, N (= 93 109) is 
the integer number of orbital cycles between the two T*, AT* (AO 
the difference between the measured T* at the N-th orbital cycle 
and its expected value, that is AT* (AO = T* 9 - (T* +Nx P orb ). 
We can assume that the correction term in Eq. [1] AP or b, is at 
most the confidence interval for P or b given by |P08t Considering 
the maximum and minimum values of P or b within its confidence 
interval, we obtain that P 01 -b lays in the interval from -2.7x 10" 11 
to +0.7 x 10~ 11 s s _1 , at lcr confidence level. 



3.2. Spin frequency secular evolution 

In the following we will derive the secular spin frequency 
derivative comparing our measurement of the spin frequ ency 
for the 2009 outburst with the spin frequency measured by P08 
analysing the XTE J1751-305 2002 outburst. Moreover, we will 
consider possible effects on the spin frequency of the two out- 
bursts occurred in the intervening seven years. To compute the 
effect on the spin frequency of the two weak outbursts between 
the 2002 and the 2009 outbursts, we assume that, during each 
outburst, the neutron star (NS) is accreting angular momentum 
L at a rate 

dL/dt = M -\JgmR^, (2) 

where M is the mass accretion rate, G is the gravitational con- 
stant, M is the NS mass, R d is the radius at which the accreting 
matter (orbiting with Keplerian speed in an accretion disc) is 
quickly removed from the disc by the interaction with the NS 
magnetic field. We assume the working hypothe sis that the mag- 
netos pheric radius R a can be expressed as (see iRappaport et al.l 
2004 and references therein) 



R., cc M~ 



(3) 



1 See http : //cxc . harvard . edu/cal/ASPECT/celmon/ 



where M is the mass accretion rate, and that the mass accre- 
tion rate is proportiona l to the X-ray flux Fx- As reported by 
iMarkwardt etall d2002l) . the light-curve of XTE J1751-305 dur- 
ing the 2002 outburst showed an exponential decay followed by 
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a sharp break after which the flux quickly dropped below de- 
tectability. Therefore we assumed to model the X-ray flux Fx of 
each outburst showed by the source with the function 



Fx = 



F X iexp(-f/r) ifO<f<f 

if t < or t > t , 



(4) 



where Fx i is the X-ray flux at the outburst peak , f is the time- 
length of the exponential decay, and r the decay time. In Table 
[3] Fx i and r are reported for each outburst. Since the 2002 out- 
burst is the only one for which a measure of the spin frequency 
derivative was possible, we will use it as a reference for the other 
outbursts. Using the above equations and hypotheses, we can de- 
rive the spin frequency derivative and, integrating over time, the 
spin frequency variation Av in an outburst of time-length t . After 
some algebraic manipulation, we find 



Av = Vi 02 T* 



1 - exp 



'Xi02 



(5) 



where t* = 7t/6, Fxi02 is the X-ray flux at the beginning of 
the 2002 outburst and 02 is the corresponding spin frequency 
derivative. The flux and the spin frequency derivative of the 2002 
outburst are F X nn = 134(7) x 10 ~ 9 erg cm" 2 s _1 (2-10 keV), 
v i02 = 0.56(12) x IP" 1 2 Hz s" 1 (seelPOsf). 
For the same outburst uV102l obtained a t = 7.1(1) days and t = 
8.5 days, whose ratio is r/t - 0.84. Adopting the same model for 
the 2009 outburst we found r = 2.4(2) days and f in the range 
1.2-2.9 days, which implies that the ratio r/f lies in the range 
0.83-2.0. 



3.2.1. Spin-down between 2002 and 2009 outbursts 

Because a significant spin frequency derivative was detected dur- 
ing the 2002 outburst (lP08l) . we considered two frequencies, at 
the beginning and at the end of the 2002 outburst, respectively. 
The frequency at the beginning of the 2002 outburst, at Tq = 
52368.653 MJD, is v 02i = 435.31799357(9) Hz. As usual the 
error is on the last digit at 1 cr level and was computed by com- 
bining in quadrature the statistical error on the spin frequency 
estimate, 4 x 10~ 8 Hz, with the systematic error induced by the 
uncertainty on the source position, 7.8 x 10~ 8 Hz. The frequency 
at the end of the outburst, which occurred about nine days after 
the beginning, when the pulsation was detected for the last time, 
is Vf 02 = 435.31799385(16), which was computed adopting the 
mean value for the spin frequency derivative given in P08j v = 
3.7(1.0) x 10~ 13 Hz s '. The error is again computed combining 
in quadrature the statistical and systematic uncertainties on vo 
(4 x 10~ 8 Hz and 7.8 x 10~ 8 Hz, respectively) with the statisti- 
cal and systematic uncertainties in the spin frequency derivative 
(7.8 x 10~ 8 Hz and 1.2 x 10~ 8 Hz, respectively), where the sys- 
tematic error on the spin frequenc y derivative induc ed by the 
uncertainty on the source position (iBurderi et al.ll2007l) is cry 
1-14 " n j_ o;„2 m i/2 tt„ „-! _ 1 £ w m-14 



sys 
-1 



: 9.6 X 10~ 14 v 3 cr" (1 + sin 2 /?) 1/2 Hz s" 1 = 1.6 X 10~ 14 Hz s 
for XTE J175 1-305. For the 2009 outburst we considered that 
V09 — V09. Averaging Eq.[5]over the outburst length f, we obtain 
an expression for the spin frequency at the beginning of the 2009 
outburst 



v; 09 = V09 - Vi 02 t* { 1 



1 - exp 



'Xi09 
7 Xi02 



(6) 



The 2009 outburst length lies in the range 1.2-2.9 days, where 
1.2 days is the time interval in which the pulsation was de- 
tected, while 2.9 days is the maximum possible length of the 



outburst (see Fig.[TJ. Assuming Fxi09 = 36 x 10 11 erg cm 2 
s , r = 2.4(2) days, and t = 1.2 Equation [6] gives v;o9 = V09 
- 0.8(3) x 10~ 8 Hz = 435.31799255(23) Hz , while considering 
t = 2.9 gives v i0 9 = v 09 - 1.7(6) x 10~ 8 Hz = 435.31799254(23) 
Hz. The difference between the two cases is, for our purpouses, 
irrelevant. 

To compute the secular spin frequency derivative we con- 
sider the total time elapsed from the end of the 2002 outburst 
to the beginning of the 2009 outburst, namely Af = 55112.0 - 
52377.653 MJD. The secular spin frequency derivative is 



(Vj 09 ~ Vf 02) 

Af 



-0.55(12) x 10 14 Hz s _1 , 



(7) 



where the 1 cr error is computed by adding in quadrature all the 
errors. Note that this value is about one order of magnitude 
higher than the secu lar spin-down rate m easured for the AMP 
SAX J 1808.4-3658 dHartman et al.ll2009l) . and of the same or- 
der of magnitude of the secular sp in-down rate measured for the 
fastest AMP IGR J0029 1+5934 <lPatrunoll2010l iHartman et all 
1201 lb IPapitto et al"1l201 lh . Even in the worst case scenario, that 
is if we computed the secular spin frequency derivative neglect- 
ing the frequency variation of all the 2002 outburst (in this case 
the total time elapsed has been computed as Af = 55112.0 - 
52368.653 MJD, namely the total time elapsed from the begin- 
ning of the 2002 outburst to the beginning of the 2009 outburst), 
we obtain v sec = -0.43(10) x 10~ 14 Hz s _1 , which is still signifi- 
cant and compatible, within the errors, with the previous value. 

3.2.2. Effect of the spin-up during 2005 and 2007 outbursts 

We now consider the fact that on 28 March 2005 and on 4 April 
2007, two other outbursts occurred. In the following we will dis- 
cuss the effect of a possible spin-up during these two outburst, 
similar to the one observed for the 2002 outburst. It should be 
noted here that there is a possibility that the 2005 outburst is 
not associated to XTE J175 1-305, since no precise position for 
the source of the outburst is available for the 2005 event. Using 
Eq. [5] we can evaluate an order of magnitude estimate of the 
spin frequency variation that may have occurred during these 
two outburst. 

For the second outburst, which occurred on 28-29 March 
2005, the peak fl ux was 19(7) X 10"" erg cm" 2 s" 1 (2-10 keV), 
and 2.4 days dGrebenev et al.ll2005HSwank et al.ll2005l) . 

For the third outburst, 4-5 April 2007, the peak flux was 
24(5) x lO" 11 erg c r rr 2 s" 1 (2-10 keV) and t 2.3 days 
(iFalanga et all 120071: iMarkwardt & Swankl 120071) . With these 
values, using Eq. [5] and adopting f/r 1, we derive Avos = 
1.5(4) x 10~ 8 Hz and Av 07 = 1.7(3) x 10~ 8 Hz, where the errors 
were computed by propagating the uncertainties on the fluxes. 
The above results are obtained in the hypothesis that t/r =a 1 . The 
result holds even if we consider, as in the case of the 2009 out- 
burst, that t = 2f. In this case the frequency changes in the two 
outbursts are Av s = 1.8(5) x 10~ 8 Hz and Av 07 = 2.1(4) x 10~ 8 
Hz. It is clear that the effect on the value of v sec , considering only 
the 2007 outburst or considering both 2005 and 2007, is negligi- 
ble. We have thus demonstrated that, within the error, the secular 
spin frequency derivative is independent on the frequency varia- 
tions caused by possible spin-up episodes during the weak 2005 
and 2007 outbursts. 



3.3. Magnetic field 

Using the derived secular spin frequency derivative, we can es- 
timate the magnetic field strength equating the rotational en- 
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Table 3. Spin frequency values for all the observed XTE J1751-305 outbursts. 



Outburst 


Vj 


Vf 


V 




T 


Av 


References 




(Hz) 


(Hz) 


(Hz) 


(lQr 11 erg cm- 2 s" 1 ) 


(days) 


(10~ 8 Hz) 




2002 


435.31799357(9) 


435.31799385(16) 




134(7) 


7.1(1) 


28(8) 


a 


2005 








19(7) 


^2.4 


1.5(4) 


b 


2007 








24(5) 


^2.3 


1.7(3) 




2009 


435.31799255(23) 




435.31799256(23) 


36(14) 


2.4(2) 


0.8(3)-1.7(6) 


d 



References. " IMarkwardt et all (|2002[). IpoI b iGrebenev et al.1 d2005h : ISwank et al.l J2005h : c iFalanga et all d2007h : 
iMarkwardt & Swank (2007); d this work JChenevez et al.l d2009h 



Notes. In this table we report, for each outburst, the spin frequency at the start of the outburst (vj, the spin frequency at the end of 
the outburst (vf), the average spin frequency (v), the flux at the peak of the outburst in the energy range 2-10 keV (Fx the flux 
decay time scale t in the hypothesis of an exponential decay, and the inferred spin frequency variation during the outburst. 



1.06 



1.04 



~ 1.02 



c 
=J 
o 
o 

"D 

CD 
N 

W 

E 




loss rate we obtain: 



0.98 



0.96 



0.94 

0.5 1 1.5 2 

Pulse phase 

Figure 5. Folded pulse profile of the three datasets in which the 
pulsation was detected. The profile is reported twice for clarity. 



ergy loss rate to the rotating magnetic dipole emission. It is 
not straightforward which expression to use to evaluate the en- 
ergy radiated by a rotating dipole. While the classical formula 
for a rotating dipole in vacuum is well known, an equivalent 
expression in presence of matter is still debated in literature. 
Gold reich & Julian] d!969l) demonstrated that NS typical mag- 
netic field strengths are strong enough to fill the magnetosphere 
with charged particles extracted from the surface, with the re- 
sult that even an aligned rotat or emits energy. W e can write the 
amount of energy radiated as (Spitkovs kvl2006l) 



-f{6) ii 1 c~ 3 co 4 , 



(8) 



where c is the speed of light, /j is the magnetic dipole moment, 
to the NS angular spin frequency, 9 the angle between the rota- 
tion and magnetic axes and f(9) a dimensionless function which 
takes into account the energy dependence from the angle 9 and 
the effects of the presence of particles in the magnetosphere. In 
vacuum f( 9) _ 2/3 sin 2 (6Q, wh ile in case of matter in the mag- 
netosphere ISrjhkovskv (2006) propose d, on the basis of MHD 
simulations, f(9) = 1 + sin 2 (#) (see Contopoulos 2007 for more 
details). Equating the irradiated energy to the rotational energy 



M = 



Ic 3 cb 



a 8.27 x 10 26 f~H9) ll 



G cm 3 , (9) 



^)4 /2 Gcm 3 . 
Assuming a 1.4M NS mass and adopting the FPS (see 



where 745 is the NS moment of inertia in units of 10 45 g cm 2 , V2 
the spin frequency in units of 100 Hz and v_i5 the spin frequency 
derivative in units of 10~ 15 Hz s _1 . Using our estimates of the 
spin frequency and its secular derivative reported in Table [TJ we 
obtain a value for the magnetic dipole strength of /i = 2.14(23) x 
10 26 / 

suming a 

| FriedmarT & Pandharipandel Il98lt Pandharipan de & RavenhalU 
1989) equation of state, we obtain a radius of R^s = 1.14 x 10 6 
cm and a moment of inertia I = 1.29 x 10 45 gr cm 2 . Under 
the assumptions above we can give an estimate of the magnetic 
field strength at the magnetic caps Bpc from the relation that 
gives the dipole magnetic field strength at the NS surface as 
a function of the magnetic dipole moment and the angle a be- 
tween the position on the surface and the magnetic dipole axis 

(a = on the magnetic cap) B - Q-i/R^ s ) Vl + 3 cos 2 ff. We find 

B PC = 3.3(4) x 10 8 f-i(6) G. Adopting f(9) = 2/3 (in line with 
what is assumed in deriving the magnetic field of radio pulsars), 
we find Bpc = 4.0(4) x 10 8 G, which is quite reasonable for this 
kin d of source. 

iHartman etaLl (1201 lb noted that the AMPs for which the 
secular spin-down was measured are suitable to be detected as 
y-ray millisecond pulsars by the Fermi Large Area Telescope, 
since several milliseco nd radio pulsars with similar characteris- 
tics were detected (see Abd o et al . 2009). The spin-down power, 
defined as E = An 2 Icoco, for th i s sour ce is £ = 0.9(2) x 10 35 /45 
erg s~'. Following Abd o et al.l d2009l) . the upper limit on the y- 
ray flux is rjE/d 2 < 2.4 X 10 33 /4s erg s" 1 kpc~ 2 , where 77 is the 
y-ray efficiency (see l Abdo et al . 2009, for 77 definition) and d is 
the XTE J175 1-305 distance lower limit, estimated by lP08l to 
be 6.3 kpc. The obse rved values for 77 lays in the range 6-100 
% (lAbdo et al.ll2009l The chance to detect XTE J175 1-305 as 
y-ray is then unlikely, although possible if the y-ray efficiency is 
high. 

For completeness, NS magneto-dipole radiation is not the 
only process to invoke for NS angular momentum loss. NS mass 
distribution can deviate fro m a perfectly sp herical distribution 
for several reasons (see, e.g. lBildstenl [l998). introducing a neu- 
tron star's mass quadmpole moment which permits emission of 
gravitational waves (G W) at a freque n cy dou ble o f the NS spin 
freque ncy. Following Hartm arTet al.l (1201 ll) and iPapitto et al.l 
(1201 ll) . it is possible to give an upper limit to the average neu- 
tron star's mass quadrupole moment Q, under the hypothesis 
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that the spin down is due only to GW emission. Using the ex- 
pression fOTjhe nettorque due to a mass quadrupole moment 
given in iThornd (Il980h and adopting the value of the spin fre- 
quency and its derivative obtained in this work, we can derive an 
upper limit to the quadrupole ellipticity (see, e.g. iFerrar l l2O10h 
Q/I < 4.6 x 10~ 9 I.y 2 (3cr confidence level), in li ne with the val- 



ues obtained for the sources SAX J1808. 4-3658 (Hartman et al 



2008) and IGR J00291 +5934 (Hartma n et al.l201 ltlPapitto et al 



201 1). It it also in agreement with the upper limits obtained for 
the millisecond radio pulsars. With these values of ellipticity and 
source distance there is no chance it can be detected by the cur- 
rent GW detectors (e.g. Abbott et al. 2010), since their sensitiv- 
ity is orders of magnitude worse. 
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